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Abstract: 

In QCD with 16| — e massless quark flavours there is an infrared fixed point with 
a s /ir = g|je, in the limit e — ► 0+. I develop the idea of Banks and Zaks to expand about 
Nf = 16±. This expansion is certainly useful for Nf = 16,15,14,..., and arguably it 
can reach the phenomenologically interesting case Nf = 2, where it suggests that a s /n 
"freezes" to a value of order of magnitude 0.4 in the infrared. 



1. In QCD with Nf massless flavours asymptotic freedom is lost if Nf exceeds 16 §, 
which is where the first coefficient of the f3 function vanishes. In my notation (5 has the 
form: 

/5(a) = ~ba 2 (l + ca + c 2 a 2 + ...), (1) 

where a = a s /ir, b = (33 — 2Nf)/6, and be = (153 — 19iVy)/12. Since b has been factored 
out, c, C2, . . . will naturally each have a simple pole at Nf = 16|. For Nf just below 16| 
there is a zero of the (3 function at a* ~ — 1/c, a small positive value. If the couplant 
a(/i) lies between and a* at some energy scale /x, then it is trapped in this range at all 
energies. One then has weak coupling at all energies [jl], Qj. 

In 1982 Banks and Zaks Q suggested an expansion in powers of 16 1 — Nf. The idea 
has since been used by White || to study the Pomeron in QCD. Grunberg Q (though in 
a context I deplore) has discovered some important features of the expansion. I argue here 
that this expansion is relevant to the real world with just two light flavours. It implies that, 
as a perturbative effect, the QCD couplant "freezes" in the infrared. This corroborates Ref. 
||, which studied R e + e - in third-order "optimized" perturbation theory ||. 

Readers who find this hard to swallow should nevertheless read on: This expansion is 
fun and it certainly offers insight into theories with 16, 15, 14, . . . flavours, whether or not 
one believes it is directly useful for Nf = 2. 

The following analysis stays entirely within the realm of perturbation theory, and I 
shall discuss nonperturbative effects only briefly at the end. 

2. The natural expansion parameter for the Banks-Zaks (BZ) expansion is 

oo = 3^-(16*-JV», (2) 

which is the limiting form of a* ~ —1/c as Nf — > 16^ from below. Because the coefficient 
8/321 is so small, ao is tiny (~ 1/80) for Nf = 16 and remains of modest size (0.36) even 
for Nf = 2. 

It is convenient to re- write perturbative coefficients, eliminating Nf in favour of ao 0]. 
The first two /^-function coefficients are: 

, 107 1 

b=—a , c = hci )0 , (3) 

with cifi = 19/4. The higher-order /3-function coefficients are renormalization-scheme 
(RS) dependent. In any 'regular' scheme their Nf dependence is such that one may write 
them as [|, f§: 

C2 = — (c2 ,-i + c 2 ,oao + C2,iao + c 2 ,2ao ) > ( 4 ) 
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and so on. (Note that a term Cij% can be assigned a degree % + j — p, and all terms in 
any formula have the same degree.) In the MS scheme one has Q: 

1 / 37117 243 34775 2 \ 
C2(MS) = ^ (-10272 + ^2- ao + W fl °)- (5) 

Consider some perturbatively calculable physical quantity in QCD. The prototypical 
example is R e + e -(Q) = 3X]o 2 (l +7^)> where 

TZ = a(l + na + r 2 a 2 + ...). (6) 

In any 'regular' RS the coefficients rj are ith-order polynomials in Nf, and hence in ao~. 

n = n,o + n,iao, (7) 

r 2 = r 2fi + r 2) ia + r 2 ^ 2 a\, (8) 

etc.. All these coefficients are RS dependent. In the MS scheme for the R e + e - case one 
has | 

n(MS(/i = Q)) = l + ^(H-8C(3))ao, (9) 



r 2 (MS(/z = Q)) 



12521 + 13C(3) 



O(ao). (10) 
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It is noteworthy that C(3) does not appear in r 10 and <^(5) does not appear in r 2 . (I have 
ignored the (^ ^) 2 term in r 2 ; its dependence depends on the electric charges assigned 
to the ficticious extra quarks.) 

3. Consider first the BZ expansion for TZ* = 7Z(Q — > 0). One first solves the fixed- 
point condition (3{a*) = for a* as a series in do- One then substitutes in Eq. @, again 
expanding in powers of do- To leading order TZ* = a* = ao, while at second order one has 

1 

a* = a [l + (c 2 _i + ci, )a + O(ag)], (11) 

and hence 

^* = ao[l + (n,o + c 2 -i + c 1)0 )a + O(ag)]. (12) 

The sum of and c 2 -i is RS invariant, as I show later, (c^o = 19/4 is invariant.) In 
the e + e~ case one obtains TZ* = ao(l + 1.22ao + ...), so the correction is relatively modest. 
At Nf = 2 the correction is about 44%. While one cannot give too much credence to the 
quantitative result, the qualitative message of leading order remains; there seems to be a 
fixed point of modest size. Figure 1 shows TZ* as a function of Nf and compares first- and 
second-order BZ results with the OPT results from Ref. M. 
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The crucial test of this interpretation will come at next order. A straightforward 
calculation yields the coefficient of the next order term in (|12|), which is 



(ci )0 + 2c 2 -i + 2ri )0 ) (ci, + c 2 -1) + ri,i + r 2 , + c 2 , + c 3 _i. (13) 

In the e + e~ case this reduces to —18.25 + 03,-1 (MS). For the expansion to be credible 
one needs 03,-1 (MS) to be in the range, say, +13 to +21. I expect 03,-1 to be found in 
the lower end of this range, thereby reducing TZ*, and bringing it into better agreement 
with the OPT results. A calculation of the 4th-order (3 function coefficient would test this 
prediction. 

Note that nth order in the BZ expansion requires n + 1 terms in the (5 function, but 
only n terms in TZ. Thus, in terms of diagrammatic information used, it is intermediate 
between nth and (n + l)th order of ordinary perturbation theory. 

The coefficients in the BZ expansion of TZ* are RS invariant. One can prove this by 
considering the RS invariants pi || and expanding them in powers of ao- From the leading 
1/ao term in p\ = p 2 + c 2 /4 = r 2 + c 2 — r\ — cv\ one sees that ri,o + c 2 ,_i is invariant, as 
claimed earlier. From the subleading part of p% one finds that ?" 2 ,o+c 2 ,o — r i o+ r i,i — ci,o r i,Q 
is invariant. Then from the leading 1/ao term in p 3 = + ±c 3 — ri(c 2 + 3r 2 — 2rf — \cr\) 
one finds that 03,-1 — 2ri,oc 2 ,_i — rf is invariant. It is then straightforward to show that 
the combination in Eq. (|l~3| ) is RS invariant. It also follows that c 2 _i + 03-1 is invariant 



4. Next, consider a BZ expansion for 7£ at a general Q. There are three preliminary 
steps. Firstly, one integrates the /3-function equation to obtain: 



b1n(u/A) = lim 



s p(x) 



(14) 



where (3{x) = f3(x)/b, and A is a constant with dimensions of mass. The constant of 
integration, C(S), must be suitably singular in the limit 5 — > 0, and I choose || 

dx 



C{8) = Pf 
Js 



s x 2 (l + ox) ' 

= - +cln(5 + cln [c| +0(5), (15) 
a 

where P (principal value) is specified because of the pole at x = — 1/c when c is negative. 
This choice amounts to a definition, in a general RS, of the A parameter. The commonly 



used A parameter [ 10 1 is defined in a less natural way, and is related by an RS-invariant, 
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but Afy-dependent factor; ln(A/A) = (c/6)ln(2 | c [ /&). The two A's become infinitely 
different in the limit Nj — > 16 i. 

Secondly, recall that for each physical quantity 1Z there is an RS invariant ||: 

Pl = bln(pL/A)-n. (16) 

The renormalization scale p cancels out because the coefficient r\ always contains a 
bln(p/Q) piece. Furthermore, the A parameter is scheme dependent in a way that exactly 
cancels the scheme dependence of the remaining part of r± [11, 01. {p\ must be regarded as 
a whole; splitting it into pieces spoils RS invariance. It cannot be written as A + Bao with 
A and B being RS invariant (TJ]. The reason is the 1/6 factor in the RS transformation 
of A: For two schemes related by a' = a(l + v\a + . . .) one has ln(A'/A) = v\/b pT[| . For 
'regular' schemes v\ is linear in Nf, but is otherwise arbitrary. The other pi invariants can 
be split into different orders in ao because they are not Q/A dependent.) One may think 
of pi as 61nQ/A e fj, where A e fj = Aexp(ri/6) is a scale specific to the particular physical 
quantity 1Z. Each A e g- can be related in an exactly known way, once the corresponding r\ 
has been calculated, to the universal A of some reference scheme (say, MS), which plays 
the role of the single free parameter of the theory. 

Thirdly, let us reconnoitre the BZ limit. Since a is trapped between and a*, it is at 
most of order ao- The first two terms in the (3 function dominate, so Eq. (|14|) gives 

TTco +OM ' (17) 

Combining the two last equations yields p\ in terms of a. The dominant terms are of order 
1/ao, so one may discard the 0{c2a) and r\ terms, which are of order unity. Since in the 
BZ limit c ~ —1/ao, ^ ~ a anci "° a * ~ a 0) tli e limiting form is: 

1 1 /TV -K\ , 10 . 
p x = - + — In — — . (18) 



61n(^/A) = - + cln 
a 



11 11* \ K 

Inverting this equation (numerically) would give 1Z as a function of p\ = b In Q/ A e g , and 
hence as a function of Q. The resulting function 1Z(Q) is RS invariant, and is universal, 
except that A e fj depends on the specific physical quantity considered. 1Z(Q) exhibits both 
asymptotic-freedom as Q — > oo and "freezing" behaviour, 1Z(Q) — > 1Z* , as Q — > 0. 

The formulation of a BZ expansion for 1Z at a general Q is not a completely unam- 
biguous matter. 1Z(Q) is not expressible as a simple power series in ao, so some thought 
is required in deciding how precisely to define the nth-order approximant. The important 
point, as with any approximation, is to reconcile and make best use of all available in- 
formation. Simply integrating the /^-function equation and then expanding in powers of 
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ao produces correction terms with (ao — TZ) denominators. Higher orders bring in ever 
more singular terms. However, these terms simply arise from an expansion of ln(7£* — TZ), 
reflecting the fact that the fixed point TZ* does not stay at ao, but is itself a series in ao- 
Therefore it is sensible to organize the expansion to reflect this. 

Thus, before performing the integration in I first re- write the l/f3(x) integrand: 

(19) 



— x 2 (l + cx + C2X 2 + . . .) x 2 (a* — x)P(x) ' 
ensuring that the pole is in the right place. Next, I express it in partial fractions: 

' 1 +Z~ ^1 ^ + H(x). (20) 



p( x ) x 2 x 7* (a* — x 
The coefficients of the first three terms are determined by the x — > and x — > a* limits. 
Hence, 7* is j*/b, where 7* is the slope of the (3 function at the fixed point: 

= -ba*(l + 2ca* + 3c 2 a* 2 + ...). (21) 



* _ dP(x) 

7 



dx 

The remainder term can be expanded as a power series, H(x) = H$ + H\x + . . .. 

In nth order of the BZ expansion one may truncate the (3 function after n + 1 terms. 
In that case H(x) = Q(x)/ P(x), where P(x) and Q(x) are polynomials of degree n — 1 
and n — 2, respectively. (For n = 1, Q(x) vanishes.) The coefficients of P(x) are of order 
unity as ao — ► 0. The coefficients of Q(x) are of order ao because of cancellations that 
make both c + I/7* and a* /j* — 1 of order ao- Thus, H(x) has coefficients of order ao- 
[Ho, for instance, is 00(04,-1 + 2e2 J -iC3 ) _i + c| _i) + 0{oq).\ In nth order (n > 4) of the 
BZ expansion one needs coefficients up to H n ^^\ for the first three orders one can drop 
H(x) altogether. 

It is now simple to perform the integration in (|14|) and use ( |l~5D to obtain p\ in terms 
of a and a*. One may then eliminate a and a* in favour of 7Z and TZ* , working to the 
appropriate order. This last step can be short circuited by noting that the final result 
must be RS invariant, and so, without loss of generality, one may choose to work in the 
RS in which TZ = a. Thus, the result in nth order of the BZ expansion can be expressed 
as: 

1 1 , A n \ , „ , H7Z i+1 



o 1 = _ + __ln 1 ) +c \ n n c \TZ) + y^ • (22) 

where TZ*^ and 7*( n ) are the nth-order approximations to TZ* and 7*, respectively. For 
small TZ (i.e., at large Q) this formula will agree with (n + l)th-order perturbation theory 
to the appropriate order in TZ and ao- For 7* a straightforward calculation gives 

f = y = ao(l + ci.oao + (c?, - c| _j - c 3 ,_i)ag + O(ag)). (23) 
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7* is the 'critical exponent' that governs the manner in which 7Z approaches 1Z* in the 



Q — > limit [13]; TZ* — TZ^ const. Q 1 . It should therefore be RS invariant. Actually, this 
statement must be qualified fll3| , [13]; 7* is invariant within the sub-class of RS's whose 
relation to the a = TZ scheme is not singular at the fixed point. The coefficients in (23) 
are the "universal" invariants (universal in that they are independent of the particular 
quantity TZ) discovered by Grunberg ||]. 

Results for TZ(Q) in first and second order are shown in Figs. 2 and 3. One could also 
obtain hypothetical third-order results by guessing a value for c& _ 1. For a value around 
15, for instance, the Q = values are similar to the first-order results, with the shape of 
the curves being more like second order. 

The first-order results use (|22| ) which, unlike the earlier form (|i~8|), retains the whole 
of c in the cln(|c| TZ) term. This seems sensible because it avoids spoiling the behaviour 
at large Q, known from ordinary perturbation theory. One must know c in order to 



obtain even the leading-order BZ result, in any case. Another point is that (22) puts 
the corrections to the coefficient of ln(l — TZ/TZ*) in the 7* denominator, rather than re- 
expanding them into the numerator. This is sensible because one knows that 7*, the slope 
of /3(a) at a*, must be positive, as is obvious from a sketch of the (5 function. 

For Nf = 16, 15, 14, . . . the theory becomes "almost scale invariant." 1Z remains con- 
stant over a huge range of Q about A e ff. This is because p\ = (107/8)ao ln(Q/A e g) 
remains close to zero. Only when Q/A e g becomes extremely small does TZ(Q) abruptly 
rise up to TZ*, while only when Q/A e g becomes very large does TZ(Q) slowly decrease, 
as required by asymptotic freedom. The constant value around A e fj is, in the BZ limit, 
TZ* /(l + x) ~ 0.787?*, where ln^ + X + 1 = 0. The region over which 1Z stays within 10% 
of this value is roughly for Q/A e g between the two extremes exp(±0.04/dg). Thus, this 
region is very extensive for clq < 0.1 (Nf > 13) and is noticeable up to ao ~ 0.2 (Nf ~ 9). 
For smaller Nf the region Q ~ A e g becomes, on the contrary, a region of rapid variation 
of 7Z(Q): asymptotic freedom sets in quickly above A e ff, while infrared "freeze-out" occurs 
just below A e ff . The cause of this different behaviour is the much larger value of the critical 
exponent 7* at smaller Nf. 

5. The above analysis has been entirely within a perturbative framework. [] Its main 
message is that the couplant "freezes" at low energies as a perturbative effect. This prop- 
erty is manifest for Nf close to 16i and the BZ expansion implies that it extends to 

1 The BZ limit might also be a useful way to explore nonperturbative effects, such as instantons, in a 
controlled weak-coupling context. 
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all lower values of Nf, albeit with diminishing accuracy. Because of "freezing" one ob- 
tains finite perturbative predictions at low energies, without resort to such notions as 
constituent-quark masses, effective gluon masses, condensates, etc.. I do not claim that 
these predictions are right — nonperturbative effects certainly do exist at low energies - 
but I do argue that these predictions are meaningful and have some predictive power. 

Let me first remark that "freezing" is not incompatible with confinement: there is 
no evidence that confinement necessarily requires the couplant to become infinite in the 



infrared. (Gribov's ideas [15], for instance, explicitly involve a "freezing" of the couplant.) 
Confinement and chiral-symmetry-breaking phenomena are associated with terms invisible 
to perturbation theory, such as exp(— 1/ba), whose Taylor expansion is + + + . . .. 
Even in the Nf = 16 theory, where the couplant is always very small, perturbation theory 
is probably not the whole story. Nonperturbative effects in that case may be extremely 
tiny but they could still be qualitatively decisive in providing confinement and chiral 
symmetry breaking. (However, if the confinement radius is enormous relative to any 
feasible experiment, then the physicists of an Nf = 16 world might well regard confinement 
as an irrelevant notion.) 

In the real world there are certainly nonperturbative effects of crucial importance at 
low energies. The perturbative predictions are smooth functions, whereas the data is 
characterized by a fine structure of hadronic thresholds and resonances. However, it is 
natural to conjecture that data and perturbative theory can be compared if some suitable 
'smearing' procedure is used to 'average out' the nonperturbative effects. This is an old 
idea (e.g. [|i|) that has appeared in many forms. I am only proposing that the game can 
be extended right down to zero energy. Indeed, such a comparison with R e + e - data proves 
to be very successful || . The hypothesis is worth pursuing because it potentially enlarges 
our predictive power. It also changes the nature of the debate about how perturbative 
and nonperturbative effects fit together: The question is no longer "How must we modify 
perturbation theory to prevent it giving nonsense for Q < A?" but instead becomes "How 
do nonperturbative effects modulate the smooth perturbative result at low energies to 
produce fine structure?" 

Quark masses were ignored in the preceding analysis. However, it should be straight- 
forward to include running masses in the usual way |17| |. The c, b, t quarks clearly decouple 
at energies below the A scale; the s quark is a borderline case; but it seems reasonable to 
treat the u, d quarks as approximately massless: Their running masses are of order 5-10 



MeV at a renormalization scale of 1 GeV IS], and although the running masses increase 
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at lower scales this rise is tempered by the freezing of a. In any case, the BZ-expansion 
results are rather insensitive to whether the effective Nf for Q < A is taken to be 3 or 2 or 
0. (It would not be appropriate to invoke constituent quark masses, which have no well- 
defined connection with the parameters of the underlying Lagrangian jl8| . They represent 
a phenomenological attempt to parametrize some of the nonperturbative effects.) 

"Freezing" has long been a popular, ad hoc hypothesis, and it has been invoked in a 
great variety of successful phenomenology. The low-energy values for a s /ir so obtained 
lie, with remarkable consistency, in the range 0.2 to 0.3. (See Ref. Q for a brief literature 
survey.) Usually "freezing" is blamed on unspecified "non-perturbative effects." The BZ 
and OPT results imply, however, that "freezing" is a phenomenon present in perturba- 
tion theory at second or third order. The BZ expansion implies that "freezing" is not 
attributable to it 2 terms; these only appear in the 0(oq) part of the r2 coefficient. It also 
implies that "freezing" is universal, regardless of spacelike/timelike Q 2 , and should occur 
for any perturbatively expandable physical quantity. 

Acknowledgements 

I am grateful to Alan White for introducing me to the concept of the 16| — Nf ex- 
pansion. I thank Jiff Chyla and Andrei Kataev for correspondence, and Frank Close, 
Poul Damgaard, Valery Khoze, Alan Mattingly, Chris Maxwell, and Graham Shore for 
discussions. 

This work was supported in part by the U.S. Department of Energy under Grant No. 
DE-FG05-92ER40717. 



8 



References 

[1] W. E. Caswell, Phys. Rev. Lett. 33 (1974) 244. 

[2] T. Banks and A. Zaks, Nucl. Phys. B 196 (1982) 189. 

[3] A. R. White, Phys. Rev. D 29 (1984) 1435; in Hadronic Matter in Collision, edited 
by J. Rafelski (World Scientific, 1989); Int. J. Mod. Phys. A 8 (1993) 4755. 

[4] G. Grunberg, Phys. Rev. D 46 (1992) 2228. 

[5] A. C. Mattingly and P. M. Stevenson, Phys. Rev. Lett. 69 (1992) 1320; Phys. Rev. 
D 49 (1994) 437; A. C. Mattingly, Rice preprint DE-FG05-92ER40717-8. 

[6] P. M. Stevenson, Phys. Rev. D23 (1981) 2916. 

[7] A. L. Kataev, in Proceedings of the QCD-90 Conference, Montpellier, 1990, ed. S. 
Narison, Nucl. Phys. Proc. Suppl. B 23 (1991) 72. 

[8] O. V. Tarasov, A. A. Vladimirov, and A. Yu. Zharkov, Phys. Lett. B 93 (1980) 429. 

[9] S. G. Gorishny, A. L. Kataev, and S. A. Larin, Phys. Lett. B 259 (1991) 144; L. R. 
Surguladze and M. A. Samuel, Phys. Rev. Lett. 66 (1991) 560. 

[10] W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta, Phys. Rev. D 18 (1978) 
3998. 

[11] W. Celmaster and R. J. Gonsalves, Phys. Rev. D 20 (1979) 1420. 

[12] W. Celmaster and P. M. Stevenson, Phys. Lett. B 125 (1983) 493. 

[13] D. J. Gross, in Methods in Field Theory, edited by R. Balian and J. Zinn-Justin 
(North-Holland, Amsterdam, 1976). 

[14] J. Chyla, Phys. Rev. D 38 (1988) 3845. 

[15] V. N. Gribov, Lund Report No. LU-TP 91-7, 1991 (unpublished); F. Close et al, 
Phys. Lett. B 319 (1993) 291. 

[16] E. D. Bloom and F. Gilman, Phys. Rev. D 4 (1971) 2901; E. C. Poggio, H. R. Quinn, 
and S. Weinberg, Phys. Rev. D 13 (1976) 1958. 

[17] H. Georgi and H. D. Politzer, Phys. Rev. D 14 (1976) 1829. 

[18] J. Gasser and H. Leutwyler, Phys. Rep. 87 (1982) 77. 



9 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402276v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402276v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402276v2 



